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ABSTRACT 

The decay rate of cosmological gravitational potential measures the deviation from Einstein-de Sitter 
universe and can put strong constraints on the nature of dark energy and gravity. Usual method to 
measure this decay rate is through the integrated Sachs-Wolfe (ISW) effect-large scale structure (LSS) 
cross correlation. However, the interpretation of the measured correlation signal is complicated by 
the galaxy bias and matter power spectrum. This could bias and/or degrade its constraints to the 
nature of dark energy and gravity. But, combining the lensing-LSS cross correlation measurements, 
the decay rate of gravitational potential can be isolated. For any given narrow redshift bin of LSS, the 
ratio of the two cross correlations directly measures [dlnD,j > /dhia]H(z)/W(x,Xs), where D$ is the 
linear growth factor of the gravitational potential, H is the Hubble constant at redshift z. W(x,Xs) 
is the lensing kernel and % and Xs are the comoving angular diameter distance to lens and source, 
respectively. This method is optimal in the sense that (1) the measured quantity is essentially free of 
systematic errors and is only limited by cosmic variance and (2) the measured quantity only depends 
on several cosmological parameters and can be predicted from first principles unambiguously. Though 
fundamentally limited by inevitably large cosmic variance associated with the ISW measurements, it 
can still put useful independent constraints on the amount of dark energy and its equation of state. It 
can also provide a powerful test of modified gravity and can distinguish the Dvali-Gabadadze-Porrati 
model from ACDM at > 2.5cr confidence level. 

Subject headings: Cosmology: the large scale structure: cosmic microwave background: gravitational 
lensing 


1. INTRODUCTION 


The decay rate of gravit ational potential is a p ow- 
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Its direct observational cons equence is the 
integ rated Sachs-Wolfe (ISW) effect l| Sachs fc Wolfel 
119671) . one class of secondary CMB temperature 
fluctuations. However, since the signal is over¬ 
whelmed by primary cosmic microwave background 
(CMB), it has to be measured b y cross correlatini 
the lar ge scale structu re (LSS ) (iCrittenden fc Turn 
1996: iSeliak fe Zaldarrialal Il999bl [), or from CMB po- 
la rization indu ced by intra-cluster electron scattering 
ifCoorav et al l 200 11. Progresses in precision CMB mea¬ 
surements and large scale galaxy surveys have enabled 
detections of the ISW-LSS cross correlation and con- 
firmed the decay of cosmological gravitational potential 
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. 2003||Afshordi et al. 
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20Q5|). Along with the flatness of the universe con- 


of the gravitational potential decay have already pro- 
vided strong evidence of the existence of dark energy 
llCorasaniti et ali 12005 1 and constraints to the dark en¬ 
ergy equation of s tate can be further improved by fu¬ 
ture observations iPogosian et ID mu. However, the 
strength of the cross correlation signal depends not only 
on the dark energy density Hde and its equation of state 
parameter w de, but also on the matter power spectrum 
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and evolving galaxy bias 1 , which can not be predicted 
from first principles. These nuisance parameters and 
modeling uncertainties could degrade the power of ISW- 
LSS cross correlation to constrain dark energy. However, 
as proposed in this paper, by the aid of gravitational 
lensing-LSS cross correlation, the evolution of gravita¬ 
tional potential can be isolated. 

The ISW effect and lensing-LSS probe the same 3D 
gravitational potential 0, with different prefactors. Since 
galaxy redshifts are observable, one can correlate ISW 
or lensing with galaxies in a narrow redshift bin. For 
such narrow redshift bin, the ISW-galaxies cross corre¬ 
lation measures <f>-5 g cross correlation, while the lensing- 
galaxies cross correlation measures (fi-5 g cross correlation. 
Here, 6 g is the over-density of galaxies. The ratio of the 
two correlations in the same redshift bin at the same 
scale then measures <//</, with prefactor which only de¬ 
pends on the geometry of the Universe, but not galaxy 
bias nor matter power spectrum. Such measurement in¬ 
volves essentially no assumptions and least amount of 
unknown cosmological parameters, so it can put robust 
constraints on cosmology and has the power to distin¬ 
guish dark energy from modified gravity. Furthermore, 
since these two cross correlation measurements are cor¬ 
related, errors in the denominator and numerator partly 
cancel. The S/N of the measured ratio is slightly better 

1 The cross correlation signal is proportional to r, the cross corre¬ 
lation coefficient between galaxy overdensity and the gravitational 
potential. In parameter fitting combining galaxy power spectrum 
measurement, r = 1 at relevant scales is implicitly assumed. But 
stochasticity could cause r to deviate from unity. So in principle, 
r should also be treated as a free parameter to be marginalized. 
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than the ISW-LSS measurement. Since cp/cp of each red- 
shift can be recovered, one obtains stronger cosmological 
constraints than that of t he ISW-lensing measurement 
llSeliak fc Zaldarriagall999ll . 

2. ISOLATING THE DECAY RATE OF GRAVITATIONAL 
POTENTIAL 

Time variation in the cosmological gravitational 
potential causes te mperature fluctuations in CMB 
ilSachs fc Wolfein 9fi7ft 

AT /' . 

-- =/[<£- ip\adx , (1) 

JCMB J 

where <p and ip are two gravitational potentials in the 
Newtonian gauge. In dark energy models, at late time, 
there is no anisotropic stress and <p = —tp. The I SW-LS S 
cross correlation power spectrum is given by (iPeebled 

[mi 

/ dk 

z = 0)—A I (k 1 l)A g (k, l) . (2) 

Here, galaxies at z — Az/2 < z < z + Az/2 
have been adopted to be tracers of the LSS. A g = 
J™ ji(kx)n{x)D g dx and A f = 2 f Q XCMB ji(!cx)-D<f,adx, 
where X 1 jX 2 >Xcmb are the comoving angular diameter- 
distance to z — Az/2, z + Az/2 and the last scattering 
surface, respectively. D g and D$ are the linear growth 
factors of galaxy over-density and gravitational potential 
<p, respectively. n g (x) is the number of galaxies per dis¬ 
tance interval. A^ s is the cross correlation power spec¬ 
trum variance between <f> and galaxy over-density <5. ji is 
the spherical Bessel function. This equation assumes lin¬ 
ear evolution in <p and galaxy over-density, which should 
be valid at large scales where almost all ISW signal comes 
from. But it does not assume the galaxy bias to be scale 
independent. Eq. [2] and El (please refer to the appendix 
for the derivation) do not require the small angle approx¬ 
imation and the Limber’s approximation and thus apply 
to all angular scales and redshift bins relevant. 

One can construct the projected gravitational po¬ 
tential $ (or its harmo n ic m od e) from C MB l ens- 
ing llSel iak fc Zaldarriagal 1999at iZaldarriaga fc Seliakl 
EMU ^^k^jnQto|j2fl02lh 21cm backgrou nd lensing 
( Zahn & Zaldarriagall2f)0rsD . cosmic shear (see iRefregied 
(2003[]_for a r ecent review) and cosmic magnification 
I Z h an g & Penl 120051) . $(n) = — f d\{(p — ip)W(x, Xs), 
where W = (1 — x/Xa)/Xi X and Xs are the comoving 
angular diameter distance to lens and source. The cross 
correlation between $ and LSS is 

r rJU 

Cq> g = 47t J Al s {k,z = 0)—A*(k,l)A g (k,l) (3) 

where A$ = 2 f 0 Xs ji(kx)W(x, Xs)D (j> dx, when <p+ip = 0. 

Eq. EH fc EH c an be further simplified by t he Limber’s ap¬ 
proxim ation llLimbeilll954 IK aise r 1998f) . lAfshordi et all 
(2004) showed that, even for a shallow survey such as 
2MASS, which has an equivalent Az ~ 0.1, the Limber’s 
approximation is still valid to several percent level for 
l > 3. We then adopt Az = 0.2 and apply the Limber’s 
approximation. Under this approximation, we have 

4^^ rj 

Cig = — A ls(-,z = o )n g D g D^a X dx (4) 
1 Jxi A 
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Fig. 1.— The accuracy of measured / = aHd In D^/d In a/ W(x), 
assuming both CMB and galaxy surveys cover the full 
sky. Fiducial cosmology has (f2 m , £2^, h, as) = 

(0.268,0.732,0.044,0.71,0.84). Transfer function is calcu¬ 
lated using the fitting formula of JEisenstein &; Hl3 119981) . 

Since D^, t)^, a, x, A| 5 (^) do not vary significantly over 
the redshift range [z — Az/2,z + Az/2] as long as that 
z > Az and n g does not vary significantly, we have 

47T 2 „ l 

Cig - -ir[ A ls(-y z = 0)n g D g D^ax]\z(X2 ~ Xi) ■ ( 5 ) 
L X 

For the same reason, we have 

C ®9 - z = 0 )n g D g D < i > Wx}\z{x2 - Xi) -(6) 

These approximations are accurate to several percent 
level for relevant l and z > 0.2. We then have 

C Ig ~ /(f)C* 9 . (7) 

Here, f(z) = [(d\nD < j > /d\na)]aH/W(x): is the key quan¬ 
tity we want to measure. Given the large cosmic variance 
in the Cj g measurement, one can safely neglect any pos¬ 
sible errors caused by this approximation. Since Ci g ~ 
f(z)C<s> g is a key relation in this paper, here we show 
another proof. The Limber’s approximation is achieved 
by taking the limit that ji(x) —> \Jtt/{21 + 1 )5d(1 + 
1/2 — x). Under this limit, for each l, only those k 
between [(l + 1/2)/x2,(( + l/2)/xi] contribute to C/ s , 
since only those k have non-vanishing A g {k,l). For 

these k, f* 1 ji(kx)D^adx ^ 0, f*° MB ji{kx)D (j> adx ^ 0 
and Aj(k,l ) ~ f* 2 ji{kx)D r j 3 adx- Similarly, we have 
A <s> - 2 f£ Mk X )W(x,Xa)D*dx for these k. Thus, we 
have Ci g ~ f(z)C$ g , from Eq. Q&E1 

/(z) is determined by the matter density, dark energy 
density and equation of state, but it does not depend 
on galaxy bias or matter power spectrum. Since Hq, 
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the Hubble constant at present, in H{z ) and W(x) can¬ 
cel each other, / does not depend on Ho, too. f{z) of 
each redshift bins can be estimated from the estimator 
/ = J2iCi g (l)wi/ J2C$gWi, where toj is the weighting 
function. For <f> reconstructed from future cosmic shear, 
cosmic magnification and CMB lensing, the fractional er¬ 
ror in the denominator is small and one can do Taylor 
expansion to estimate the error in /. For this limit, we 
have 


A/ 2 = E i Agf g wf J2i AClgwf 

P E l C IgWl } 2 [E i C *g w l? 

Here, 


2 E; (AC/ g A$ g )u; 2 

E; ClgWl E l C^gWl 
( 8 ) 


A C 2 g {l) = 


A Cl g (l) = 


CcmbCq 




{21 4 

C*Cg 


l)/sky 

hC 2 s 


{21 + 1) ,/sky 


(9) 


{AC Ig AC^ g ) = 


C^jCg + ClgC&g 


{21 + l)/sky 

are the statistical errors of C/ s , (7$ ff and cross correla¬ 


tion between them, respectively. Ccmb, C g and (7$ are 
the sums of corresponding signals and associated con¬ 
taminations. For <f> reconstructed from CMB lensing of 
a very low noise CMB experiment, (7$ ~ (7$ at l < 200 
llHu fe Okamoto! 1200 2*). For $ reconstructed from the 
cosmic magnification and cosmic shear of future radio 
and optical surveys, this holds over even larger l range. 
Since m ost ISW-LSS cr oss correlation signal comes from 
l < 100 tlAfshordill2004D and 21cm emitting galaxies have 
high surface density, we neglect shot noise term 2 and ap¬ 
proximate C g = C g . Since £7 CMB /(7 ISW i ; errors in 
Gjg will dominate over errors in C^ g and cross correla¬ 
tions. Furthermore, the last two terms in Eq. 0 partly 
cancel. So it is a good approximation to neglect the last 
two terms. Under this simplification, we obtain the min¬ 
imum variance estimator wi = Pzzrjis. and the minimum 

‘ XL 'Ig\ L ) 

variance is 


A/ 2 Cj g \ (2 1 + l)/sky<7 ISW \ 

P ~ y , AC'jgJ r 2 CC mb J 

( 10 ) 

where <7 ISW is the ISW power spectrum of the corre¬ 
sponding redshift bin and r is the cross correlation co¬ 
efficient between the ISW effect of the corresponding 
redshift bin and galaxies. For narrow redshift bins we 
choose, r is very close to unity. 

We choose the fiducial cosmology with best fit WMAP 
parameters = 0.268, fl\ = 1 — fib = 0.044, 


2 For the estimations of LSS clustering signal and shot noise, 
biggest uncertainties are (1) HI (neutral hydrogen) mass function 
at high z, (2) 21cm emitting galaxy bias and (3) specifications 
of 21cm experiments. If one adopts HI mass functions calibrated 
against observations of damped Lyman-a systems and Lyman limit 
systems, SKA can detect ~ 10 9 galaxies at z ~ 3 in five years 
across the whole sky, for a field o f view 10 deg 2 at ^ 3 00 Mhz (for 
details of the calculation, see, e.g. IZhang &; Penl J200(3H . For SKA, 
detection threshold of HI mass at z = 3 is ~ HFMq, so detected 
galaxies are likely having biases bigger than one. Then, one can 
neglect the shot noise term with respect to C g 
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Fig. 2. — ^de-^de contour. Since most signals come from 
Zif. ~ 0.5, d In Dj, /din a at z* ~ 0.5 should be close to the fiducial 
cosmology in order to get a good fit. This requires that the dark 
energy density Qde( 1 + £*) 3 ( 1+ ™ DE ) ~ Q^ ducial . We plot the line 
with z* = 0.5. 


h = 0.71, as = 0.84 and primordial power index 
n = 1. Transf e r fun ction is calculated according to 
lEisenstein fc Hd 1 1998 b We use this fiducial cosmology 
to estimate the S/N of / to be measured from future 
observations. We then use this fiducial S/N to forecast 
dark energy and gravity constraints. We choose 21cm 
emitting galaxies as tracers of LSS. 21cm emission is 
less affected by dust and 21cm surveys such as SKA can 
cover the whole sky. They also have the advantage to 
go deep into z > 3 and allow the measurement of / at 
z ~ 3. We choose the CMB re ference experiment dis¬ 
cussed in iHu fc Okamoto ( 120021) to produce a full sky 
lensing map. Galaxy bin size is chosen to be A;? = 0.2. 
For each galaxy redshift bin, one obtains a measure of 
f(z) at that redshift. The result is shown in Fig. 0 
For the fiducial model, the decay of </> can be detected 
at 2 < 2. Since the physical separation of two redshift 
bins are much larger than the galaxy correlation length, 
LSS of any two redshift bins are uncorrelated, / of each 
bins are close to be independent 3 . We will do a reduced 
X 2 analysis to estimate the dark energy constraints and 
gravity constraints. The reduced y 2 is defined as 


a* 2 = E 

i 


if Pi) - /f d ) 2 


A/, 


fid.2 


( 11 ) 


where f{zi) and / fid are the predicted / of a given dark 


3 Because the ISW effect is projected over the redshift range from 
zero to ~ 1100, there does exist a correlation between Ci g (l) of dif¬ 
ferent redshift bins. The same is true for C$ g . However, statistical 
errors of / in each redshift bins are dominated by primary CMB. 
So the cross correlation coefficients between error bars of different 
redshift bins are close to zero. It is safe to neglect such correlations 
and take / measurements of each redshift bins as uncorrelated. 
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energy model or gravity model and of the fiducial model, 
respectively. 


3. CONSTRAINING DARK ENERGY 


The quantity / is very sensitive to the amount of dark 
energy. To visualize it, we plot those of f2 m = 0.2, 0.4 flat 
ACDM in Fig. |T| One can see that these two cosmologies 
can be distinguished from the fiducial cosmology at > 3cr 
confidence level (C.L.). Given a flat ACDM prior, Da can 
be constrained to 0.73lo!o7 at 2er C.L.. 

We want to constrain both Ode and icde simultane¬ 
ously. At super-horizon and near horizon scales, dark en¬ 
ergy fluctuation can contribute ~ 10% to the decay of </> 
dBean &; Dordl2004 iHu k, Scrantonl feOOd') . To calculate 
the dark energy fluctuation, one extra parameter, the 
sound speed c s , is required. Fortunately, the ISW signal 
(weighted by noise) mainly comes from sub-horizon scale. 
For the purpose of this paper, we can neglect the dark 
energy fluctuation in the parameter estimation. With 
this simplification, D$ is given by 


D" I D ( 5 I H ) I Drk (3 [ H a 3 ) - Q (12) 

* + * [ a + H )+ a 2 H 2 a 3 H 2 ’ ° ’ { 2) 

where = d/da and = d 2 /da 2 . We numerically solve 
this equation with the initial condition that when a —■> 0, 
—>const. and D^ —> 0. The results is shown in Fig. 
[3 The constraints are not very impressive, comparing to 
other methods. But since the theoretical prediction only 
involves general relativity and parameter fitting involves 
the least amount of free parameters, the constraints are 
less affected by model uncertainties and the existence of 
dark energy can be confirmed at high confidence level. 


4. DISTINGUISHING DARK ENERGY FROM MODIFIED 
GRAVITY 


Precision cosmology provides crucial tests of general 
relativity at cosmological scales. Modifications to general 
relativity change not only the expansion history of the 
Universe, but also the LSS. For modified gravities which 
reproduc e the expansion h i story of dark energy mod ¬ 
els (e.g. iDvali et al.l (l200Gf) : iCarroll et"afl {2004 l2005h: 
iMena et alJ 11200,4) 1. the structure growth is in general 
different. Thus, the measured / from our methods can 
provide crucial tests for such models. In general, modi¬ 
fications to general relativity causes to be scale de¬ 
pendent. In this case, one has to choose much narrower 
l bin to isolate f(z) of each scale. Since the relative 
error in the denominator can approach unity for nar¬ 
row l bin, more sophisticated estimator is required to 
isolate /. For simplicity, in this paper we o nly d iscuss 
the D vali-Gabadadze-Porrati (DGP) model (|Dva,li et a 1.1 
2000), which preserves the scale independ ent 
llLue et alll2004bi: iKovama fc MaartensfeOO flf. 

For a flat DGP model, the expansion history is given 
by H 2 = ^ + Hgfl m a~ 3 , where r 0 = l/[#o(l - ^m)]- 
The two Newtonian potential </> and i/ no longer follow 
the relation </>+?/> = 0. The evolution oi cfy—ip is governed 

by 


D " +D ‘ ,5 H 

H )+ n2 (.<*+ H 2 a 3 H 2 G 


H 


(13) 


w here G R ff/G = 1+1/3/3 and 0 = l+2mH 2 /(l—2 r n H) < 
0 llLue et ah ' 2004b ; IKovama fc Uaaricns 200' ):. This 



Fig. 3. — Distinguishing DGP from the fiducial ACDM cosmol¬ 
ogy. Since at high redshifts gravitational potential decays much 
faster in DGP than in ACDM cosmology, DGP can be distinguished 
from ACDM at > 2.5cr confidence level. 


equation holds at sub-horizon scales with k aH and 
k 1/ro. At z > 2, l < 10 modes do not satisfy these 
conditions. But at these redshifts, most signal comes 
from l > 10 and throwing away l < 10 modes does not 
degrade the / measurement significantly. For simplicity, 
we neglect this complexity. 

We find that flat DGP model can be distinguished from 
the fiducial ACDM cosmology at > 2.5cr G.L.. Fig. 1 
demonstrates this point. When a — > 0, the expansion 
history of a DGP model resembles a u>de = —1/2 dark 
energy model (or ACDM with curvature). But gravi¬ 
tational potential in DGP decays faster, due to extra 
suppression caused by smaller effective Newton’s con¬ 
stant G e ff. Comparing to the fiducial ACDM model, a 
low D m DGP model causes too fast gravitational po¬ 
tential decay at high redshifts while a high D m DGP 
model causes too slow gravitational potential decay at 
low redshifts. So it is hard for DGP model to repro¬ 
duce / and its redshift dependence in the fiducial ACDM 
cosmology. One can infer similar conclusion From Fig. 
3, where icde > —0.6 dark energy model is excluded 
at > 90% confidence level. Since gravitational poten¬ 
tial in DGP decays faster than a icde = —1/2 dark 
energy model, DGP should be distinguished from the 
fiducial model with higher C.L.. This can be better un¬ 
derstood from the asymptotic behavior of D^. When 

a —> 0, din D^-^/d In a —► — [11/16] [(1 — fi ra )/Slm 2 ]o 3 / 2 
for DGP and dInD^-^/dlna —> —[3(1 — icde)/(6 — 
5w’de)][(1 — Hm)/r2m]a -3u,DE for dark energy. At high 
redshifts, gravitational potential decays faster in DGP 
than in dark energy models with icde < —1/2 and much 
faster than in ACDM cosmology. It is this feature that 
allows a clear discrimination between DGP and many 
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dark energy models using the ISW effect. 

5. DISCUSSION 

We have shown how to isolate the decay rate of grav¬ 
itational potential from ISW-LSS and lensing-LSS cross 
correlation measurements. The measured decay rate, 
with prefactors only depend on geometry of the uni¬ 
verse can put robust constraints on dark energy and 
the nature of gravity, free of many theoretical uncer¬ 
tainties. The accuracy of such measurement is domi¬ 
nated by statistical fluctuations dominated by primacy 
CMB. This allows us to safely neglect several possible 
svstematics. (T) Dark energy fluctuations. As shown in 
iBean fc Pore ( 2004), such fluctuations are at most 10% 
those of dark matter fluctuation. For sound speed close 
to unity, these fluctuations effectively vanishes at most 
scales we are interested and the overall effect can be ne¬ 
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